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ABSTRACT 



Sound fields in wedge-shaped ocean layers, modeling con- 
ditions on the continental shelf, have been studied at the 
Naval Postgraduate School in the last few years using the 
method of images. These studies are carried further in the 
present work. The method is implemented in different en- 
vironmental conditions. This thesis examines the influence 
of several parameters on the sound field for downslope pro- 
pagation in a wedge-shaped fluid of speed of sound C 2 over- 

lying a slow bottom of speed of sound cj.. On the basis of 

gualitative and semi-guantitative analysis of the behavior 
of the pressure-depth profile for various geometrical and 
physical parameters, we can conclude that: 

1. A defined distance, the "characteristic distance" 

X Q = ir/( 2k2 sin 0 Q tan 6), where cos 0 O = c^/C 2 , 

k 2 = w/c 2 , and 0 is the vertex angle of the wedge, 

has physical meaning as a useful scaling distance. 

2. The distance of the source from the apex, in terms 
of the X Q , plays a major role in determining the 
downslope sound field. 



4 



BUT'LEY ' 7 TT^r VP~ 

XiAV/,1 T j ■ F,*_ <" ■ ,r ^-riT, 

lion, l 7 : • -• ... .^>.o-600S 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

A. SOUND PROPAGATION IN THE WEDGE 11 

B. THE METHOD OF IMAGES 13 

C. A COMPUTER PROGRAM DSLOW 14 

II. THEORY 18 

A. GENERAL VIEW OF PRESSURE DISTRIBUTION IN 

THE WEDGE DOWNSLOPE 19 

B. SOUND PRESSURE AT A POINT IN THE WEDGE 



DOWNSLOPE PREDICTED BY THE METHOD OF IMAGES ... 28 



III. DSLOW PROGRAM IMPLEMENTATION 28 

A. PROGRAM FEATURES 28 

B. NORMALIZATION 29 

C. PROCEDURE 31 

D. PROGRAM IMPROVEMENT 31 

IV. DISCUSSIONS 45 

A. GRAPHS OUTPUT 45 

B. GRAPHS CLASSIFICATION 46 

1 . Type 1 Curves 46 

2. Type 2 Curves 46 

3. Type 3 Curves 48 

C. TRANSITION POINT 48 



5 



D. PARAMETER VARIATIONS 53 

1. Variations of 6 53 

2. Variations of y 53 

3. Variations of c^/c 2 and p^/p 2 54 

V. CONCLUSIONS AND RECOMMENDATIONS 55 

A. CONCLUSIONS 55 

B. RECOMMENDATIONS 56 

APPENDIX A: DSLOW ALGORITHM 57 

APPENDIX B: DSLOW PROGRAM 60 

APPENDIX C: NUMERICAL RESULTS OF DSLOW 65 

APPENDIX D: TABLES 71 

APPENDIX E: GRAPHS OF Rj_ VERSUS R 2 AT THE FIRST 

TRANSITION POINT 74 

LIST OF REFERENCES 82 

INITIAL DISTRIBUTION LIST 84 



6 



LIST OF TABLES 



1 . 


Receiver Distance at 
for Constant p^/ P 2 = 


the 

0.80 


2 . 


Received Distance at 
for Constant p^/ p 2 = 


the 

0.80 


3 . 


Receiver Distance at 
for Constant p^/ P 2 = 


the 

0.90 



First Transition Point, 

c l/ c 2 = 1.10 71 

First Transition Point, 

c l/ c 2 = 1.20 72 

First Transition Point, 

c l/ c 2 = 1.10 73 



7 



LIST OF FIGURES 



1.1 Geometry of the wedge 17 

2.1 Ray tracing in the wedge downslope 20 

2.2 Geometry of a wedge by the method of images 22 

2.3 Geometric development of reflection angles 24 

2.4 Geometry of symmetric images 26 

3.1 Pressure amplitude normalization 30 

3.2 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 33 

3.3 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 34 

3.4 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 35 

3.5 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R} varied 36 

3.6 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, varied 37 

3.7 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 38 

3.8 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, P± varied 39 

3.9 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R]_ varied 40 

3.10 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 41 

3.11 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, varied 42 

3.12 Graphs of receiver angle 6 versus pressure 

amplitude with R 2 fixed, R^ varied 43 



8 



3.13 Graphs of receiver angle 5 versus pressure 

anplitude with R 2 fixed, R ^ varied 44 

4.1 The plots where there are two transition 
points, R 2 = 4.6 is the first and R 2 = 6.4 

is the second 47 

4.2 Tvne 1 curves, indicating a pressure anplitude 

nearly linear with depth SO 

4.3 Tyne 2 curves, indicating a well-colligated 

sound field as the source away fron the apex ...... 51 

4.4 Type 3 curves, indicating the presence of 

reflection and diffraction near the botton 52 

E.l R^ vs P .2 at the first trans. points, for 

3 = 6°, P 1 /P 2 = 0.80, c^/c 2 = 1.10 74 

E.2 Ri vs ?2 at the first trans. points, for 

3 = 10°, p 1 /p 2 = 0.80, c 1 /c 2 = 1.10 75 

E.3 vs R 2 at the first trans. points, for 

3 = 15°, p! /P 2 ~ 0.80, c^/c2 = 1.10 76 

E.4 p-^ vs R 2 at the first trans. points, for 

3 = 6 °, p ± / p 2 ~ 0.80, C 1 /C 2 = 1.20 77 

E.5 R^ vs P 2 at the first trans. points, for 

3 = 10°, P1 /p 2 = 0.80, c 1 /c 2 = 1.10 78 

E.6 Rj vs R 2 at the first trans. noints, for 

3 = 15°, Pl /p 2 = 0.80, c 1 /c 2 = 1.10 79 

E.7 Pj vs P ~2 at the first trans. points, for 

3 = 6°, pi/p 2 = 0.90, c^/c2 = 1.10 80 

E.3 R} vs P -2 at the first trans. points, for 

3 = 10°, p 1 /p 2 = 0.90, c 1 /c ? = 1.10 81 



9 



ACKNOWLEDGEMENTS 



I would like to acknowledge Professor Alan B. Coppens 
who assigned this interesting topic. My thanks to Commander 
Chil Ki Baek who assisted Professor Alan B. Coppens thus 
making this job much easier. I thank my wife Tuti and my 
children, Henry, Anna, and Michael, for their patience and 
support. I thank those people in the Physics Department of 
the Naval Postgraduatre School who encouraged me to 
persevere . 

Finally, I hope that in the future this work will be 
useful to others by expanding their knowledge in this field. 



10 



INTRODUCTION 



I . 

A. SOUND PROPAGATION IN SHALLOW WATER CHANNEL 

Experimental investigations of sound propagation in 
shallow water channels have been done by several investi- 
gators. Shallow water propagation is of interest because of 
the applications to coastal defense. These investigations 
are expensive and time consuming. The use of a computer 
model should provide a relatively inexpensive alternative to 
observation . 

One of the techniques uses normal mode theory. The 
normal node theory, which was introduced and developed by C. 
L. Perkeris [Ref. 1], gave the exact solution in water of 
constant depth. Further development of normal mode theory 
was made by L. Brekovskikh [Ref. 2], who initiated pressure 
as an integral involving Bessel functions and solution of 
the normal mode equation. Another theoretical approach to 
sound propagation in a horizontally stratified ocean of con- 
stant depth is given by the method of multiple scattering 
(Ref. 1], With this method, all the previous theories can 
be simplified by conversion into an asymptotic form which is 
valid when the acoustic wavelength is small compared to the 
distance over which the sound speed varies appreciably. 
These theories agree with the laboratory experiments. 
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For a water channel with a small bottom slope, the sound 
field may be expressed approximately in terms of adiabatic 
normal modes. To facilitate prediction, R.D. Groves, Anton 
Nagl, H. Uberall, and G. L. Zauer [Ref. 3] modeled a wedge- 
shaped isovelocity ocean with a linearly-sloping, perfectly- 
rigid ocean floor using adiabatic normal modes. For a pene- 
trable bottom the normal mode description fails when modes 
propagating upslope encounter the "critical depth" (H c ), 
defined as the depth where the associate mode changes from 
fully trapped within the water channel to radiating energy 
into the bottom (cut off) [Ref. 4-6]. The parabolic egua- 
tion can be used to explain the mechanism of sound energy 
radiation into the bottom [Ref. 7,8] . Such an eauation was 
studied by F. B. Jensen and W. A. Kuperman [Ref. 9], with 
predictions that satisfactorily agreed with the experimental 
results for small ray angles. With some restrictions, 
normal mode theory is applicable for sound propagation in 
the wedge-shaped fluid with a fast bottom. The parabolic 
eguation is good for fast and slow bottom, but with the 
restriction that horizontal ray angles must be less than 
20 ° . 

Another technigue introduced to predict the propagation 
of sound in the wedge is the method of images. This method 
was derived from the simplest case; a monofrequency point 
source in a homogeneous ocean with parallel boundaries. The 
total pressure is the sum of an infinite number of spherical 
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waves from an infinite set of images. The restriction of 
this method is that it does not generalize to the case of 
inhomogeneous media or non-planar boundaries. In this work, 
this method will be studied. 



B. THE METHOD OF IMAGER 

In 1978, Coppens, Sanders, Ioannou, and Kawamura [Pef. 
10] , predicted the pressure amplitude and phase of the sound 
field along the bottom of a wedge-shaped fluid layer of 
density > and speed of sound Cj_, overlying a fast fluid 
bottom of density p 2 > and speed of sound c 2 > by applying 
the method of images in a computer program implementation. 
In 1934, Baek [Ref. 11], and LeSesne [Ref. 12], implemented 
further improvements. Baek's computer program, WEDGE, and 
LeSesne's computer prooram XSLOPE were validated for several 
cases. WEDGE was developed for two-dimensional upslope pro- 
pagation (the source and received are in the same vertical 
plane perpendicular to the shore line, and the receiver is 
closer to the apex than the source (Figure 1.1a)) and down- 
slope propagation (the source is closer to the apex than the 
receiver (Figure 1.1b)). XSLOPE was developed for upslope, 
downslope, or cross-slope propagation (the source, receiver, 
and apex, are not necessarily in the same plane perpendicu- 
lar to the shore line (Figure 1.1c)). In both programs, 
Baek and LaSesne assume that the fluid in the wedge and 
fluid in the bottom have constant densities, that 
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the speed o£ sound is constant, and that the interface be- 
tween the fluids and the surface is smooth. 

In both WEDGE and XSLOPE, all distances are scaled in 
units of the "dump distance." A dump distance X, as stated 
in Reference 10, is the distance from the apex measured 
alonq the interface at which the lowest mode attains cutoff. 
If the wedge angle is 8 (Figure 1.1c), then 

A ki sin0 c tan 6 ' 1 * 1 1 

6 C = arccos ( c 1 /c 2 ) (1.2) 

where k^ is the wave number in the wedge and 0 C is the crit- 
ical grazing angle for reflection of sound from the bottom. 
For 3 < < 1 

X = H tan 8 (1.3) 

This scaling distance negates the necessity of specifying 
f reguency . 

C. COMPUTER PROGRAM DSLOW 

At the start of the work reported in this thesis, a com- 
puter program was obtained (Ref. 13], which is an extension 
of the WEDGE and XSLOPE for downslope configuration with a 
slow bottom. The computer model, DSLOW, developed to run on 
a desktop computer (Wang 2000), uses the method of images to 
predict the pressure amplitude and phase anywhere 
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within the wedge fluid overlying a slow bottom in a cross- 
slope configuration. A geometrical picture of this configu- 
ration is shown in Figure 1.1c. 

Mathematically, the model used in WEDGE and XSLOPE is 
applicable in any condition. But consideration must be 
given for making it work for a slow bottom. In the case of 
a fast bottom, the dump distance has a physical meaning. 
The dump distance is expressed as a function of the critical 
angle. The critical angle is equal to arccos (ci/c 2 ). In 
the case of slow bottom, cj_/c 2 is greater than 1, thus 
arccos (c^/c 2 ) is invalid; therefore, so is the dump 
distance. To facilitate the scaling factor, a "characteris- 
tic distance" or "scaling distance" is introduced. We need 
the scaling distance because, with this distance, our model 
will be independent of frequency as in the fast bottom case. 
There is also the hope that the use of a scaling distance 
will allow systematic observation of the pressure field. 
This scale distance X Q is the distance measured along the 
interface from the apex to the point where the lowest mode 
would attain cutoff if the fluids in the wedge and in the 
bottom were to be interchanged. The characteristic distance 
is defined by the following equation; 



X 



o 



tt/2 

k 2 sin0 o tan B 



( 1 . 4 ) 
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where 9o = arccos ( c 2 /c ^ ) and K 2 = u/c 2 is the wave number 
in the bottom. 

The following terms will be used throughout (see Figure 

1 . 1 ): 

3 = wedoe ancle 

= distance of the source from the apex in units of 
*o 

R 2 = distance of the receiver from the apex in units of 

'<o 

Y = angle of elevation of the source above the bottom 

5 = ancle of elevation of the receiver above the 

bottom 

Y 0 = distance between the projection of the source and 
receiver on the shore line, scaled by X Q 

p^/p 2 is the ratio between the density of the fluid in 
the wedae (p^) and the density of the fluid in the 

bottom ( p 2 ) 

c l/c 2 is the ratio between the speed of sound in the 
wedge (c^) and the speed of sound in the bottom 
(c 2 ). A fast bottom occurs when c 2 > c^; a slow 
bottom occurs when c 2 < c-^ 

The purpose of this research is the following: 

1. To transfer, test, and evaluate DSLOW program on the 
IBM 3300; 

2. To obtain numerical and graphical output for a 
number of cases; and 

3. To attempt to develop plausible explanations for anv 
significant features observed. 
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c. cross-slope configuration 



Figure 1.1 Geometry of the wedge 



II. THEORY 



A. GENERAL VIEW OF A WEDGE PRESSURE DISTRIBUTION IN THE 

DOWNSLOPE CONFIGURATION 

A general picture of the sound energy propagation within 
the wedge in downslope direction is given in Figure 2.1. If 
a sound source is placed at point S, ray 1 will reach the 
surface at point P with an incident angle a with respect to 
the normal to the surface at this point. This ray is re- 
flected by the surface at the same angle but with the phase 
180° different. (On the surface, sound pressure is zero 
everywhere.) The reflected ray reaches the bottom with an 
incident angle 8 + a . At great enough distance, ray 1 
never reaches the bottom again. This ray does not contri- 
bute to a sound pressure field at the bottom. The pressure 
at the bottom should be very small according to the ray 
theory argument. 

Using these ray-tracing methods, an estimated profile of 
the pressure amplitude versus the receiver depth can be 
made. When the source and the receiver are placed near the 
apex, the pressure amplitude is zero at the surface, a maxi- 
mum somewhere within the wedge, and greater than zero at the 
bottom. In the case where the source is at a far distance, 
the pressure amplitude is egual to zero at the surface, a 
maximum somewhere within the wedge, and zero at the bottom. 
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Ray tracing will only give a rough approximation, not an 
exact solution, but ray tracing may be used as a guide. The 
method of images calculates the exact pressure amplitude at 
each point within the wedge subject only to the assumption 
inherent in using the plane-wave Rayleigh reflection 
coefficients . 

B. SOUND PRESSURE AT A POINT IN THE WEDGE DOWNSLOPE 

PREDICTED BY THE METHOD OF IMAGES 

Let the source be a scaled distance R 3 from the apex and 
at an angle of y measured from the bottom of the wedge. Let 
the receiver be a scaled distance R 2 from the apex and at an 
angle 6 measured from the bottom. 

Using Figure 2.2, let the upper half family of images be 
n = 1,2,3,4,... and the lower half family be n' = 
1,2, 3, 4.... Calculating the field resulting from source and 
images proceeds along the lines developed in [Ref. 14] . If 
<j> n is the angle formed at the apex between the n t ^ image of 
the source and the receiver, then 

4>1 = 26 - 6 - y 
<j>2= 26“ 6 + y 
$ 3 = 46 - <5 - Y 
$ 4 = 46 - 5 + Y 
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Figure 2.1 Ray tracing in the wedge downslope 
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or more generally 

<ji n = (n+l)S-6-Y for n odd 
<j> n = n S - 5 + y for n even 

Which can be reduced to: 

4 > n = (n+(l/2) [1— (— 1 ) n ] } 6 + ( -1 ) n y - 5 (2.1) 

or 

4> n = 2 INT [ ° + — ]B + ( -1 ) n y - 6 (2.2) 

where INT[ ] denotes the largest integer which is equal to, 
or smaller than the argument. Using the same method for the 
member n' of the lower family of images we obtain: 

4> n « = (n+ ( 1/2 ) [l-(-l) n ] }$ + ( - 1 ) n y + 5 (2.3) 

or 

4> n t = 2 INT [ n j - ] S + ( -1 ) n y + 6 (2.4) 

Using the geometry of Figures 2.2 and 2.4, the distance 
between the n* 1 * 1 and n'^ images to the receiver is 
respectively 

r n = / + R 2 2 - 2R;[R 2 cos <}> n (2.5) 

and 

r n 1 = /R + R 2 2 - 2R 1 R 2 c °s 4> n (2.6) 
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Figure 2.2 Geometry of a wedge by the method of images 
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The angles 


9 no 


and 


0 n i o for the 


n tR and n 


» t h 


images 


respectively are 




















r sin <j> n 








9 no 


arc 


tan 


R 2 //p l “ cos ‘t'n 


] 




(2.7) 


and 
















arc 


tan 


f ? in4> n ' 

L 7K r 


] 




(2.8) 


o 






R 2 /R 1 - cos <f> 


n’ 






Define 0 nn , 


and 


9 n ' : 


n as the angles 


of incidence 


for the 


m^ bounces from 


the 


bottom for the n and 


n ' 


image 


respectively; m 


= 1 


,2,3 


(The 0 th 


bounce is 


the last 


one before reaching 


the 


receiver . ) The 


geometry 


of 


Figures 


2 . 2 and 2.3 aive 


9 nm 


as 


follows : 








e 21 = 


8 20 


- 26 


- 6 








9 32 = 


CD 

U) 

o 


- 46 


- 6 








0 4 1 = 


9 40 


- 26 


- 6 








e 52 = 


9 50 


- 46 


- 6 









The general expression is 

®nm = e no ~ 2mg - 6 
Using the same method 

e n’m = 9 n'o " 2m 6 + 6 

The maximum number of bottom bounces of the n*-*" 1 and n ,ta 
image is 

m p , ax = M = I NT [ <j> n /2 6 ] = INT [<|> n ./26] 
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Figure 2.3 Geometric development of reflection angles 
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The maximum number of images is 



n max = N = INT U/0] 



The reflection coefficients for the n*-* 1 and n'^ 
for a plane wave are: 



and 



R ( 9nm ^ 



Pl c l 
P2 U 2 

+ 

P 2 c 2 



- * 



nn 



^nm 



^ ( ®n 1 m ^ 



Pl c l 

p 2 ^ 2 ^ ^ 



p i c i. + 

P2 C 2 



♦ n' 



m 



where 



tnn 



/ 1 - (C|/C ? ) 2 COS 2 0~^~ 

sin 0 nm 



and 



^n ' m 



/ 1 - ( c i /c •p ) ^cos ^ 8 n 'm 
sin 6 n . m 



The contribution from the upper family of images 



N 

Pu = l 



M 



n= 1 rn 



exp(-jkr n ) (-1 ) INT [ n+ l )/2] n R 



m= 0 



nm 



images 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

is 

(2.13) 
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Figure 2.4 Geometry of symmetric images 
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and for lower family of images is 



Pu 



N 

- I ' 

n= 1 



exp(-jkr n » )(-l) IM T[n+l)/2] 



M 

n R n ' m 
m= 0 



( 2 . 14 ) 



The total complex pressure is 



P(x) = Pu + PI 



( 2 . 15 ) 
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Ill . 



DSLOW PROGRAM IMPLEMENTATION 



A. PROGRAMS FEATURES 

Since the mainframe graphics computer was available, the 
DISSPLA graphical program was used. The only programming 
language compatible with DISSPLA is FORTRAN. The numerical 
and graphical output is provided by this program. To give 
the pressure amplitude versus received angle graphs, two- 
dimensional plotting is used. 

The program DSLOW is run by placing the point source 
anywhere in the wedge and then placing the receiver at a 
distance downslope from the source. The receiver position 
was varied from zero degrees at the bottom to 8 at the 
surface. High resolution plotting was achieved by dividing 
the y-axis (received angle) into two regions. The first 
region covers the receiver angles from zero to 1/5 of the 
wedge angle. In this region A6 is equal to 8/100. The 
second region covers the remaining wedge angle with A6 equal 
to 8/10. This method provides 29 predictions of the pres- 
sure amplitude. Another method of plotting carried out was 
in the region of 6 > 6/2, A6 = 6/10, and in the region of 6 
< 6/2, A 6 = 6/100. This method provides 54 points to be 
plotted . 
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B. NORMALIZATION 



The main goal of this research was to investigate the 
profile of the pressure amplitude as a function of a number 
of variables. An example of the numerical values of the 
pressure amplitude, the normalized pressure amplitude, and 
the phase at each receiver position is displayed in Appendix 
C. The sound pressure becomes smaller as the receiver is 
moved away from the source. If the pressure amplitude were 
olotted directly, it would be difficult to compare the 
curves at near distances to the curves at far distances 
since at the near distances the pressure amplitude is much 
greater than the pressure amplitude in far distance. Thus, 
a normalized pressure amplitude is needed. The normalized 
pressure is obtained as follows: (see Figure 3.1) 

We know that the sound pressure at the surface is zero 
and that the sound pressure is a small number greater than 
zero at a point near the surface. The first non-zero value 
of pressure P]_ is at the receiver angle, 6-^ = 90/10. We use 
this first calculated non-zero pressure amplitude as the 
normalization unit. The normalized pressure is 

PN = P( 6)/(P 1 ) (3.1) 
where P(5) is the pressure at any point within the wedge. 
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Figure 3.1 Pressure amplitude normalization 
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C. PROCEDURE 

Figures 3.2 through 3.7 represent the results when the 
receiver distance and source angle are fixed and the source 
distance and receiver angle are varied. Figures 3.8 through 
3.13 represent the results when the source distance and 
angle are fixed and the received distance and angle are 
varied. These cases will be the foundation of our 
subsequent discussions. 

The solid lines indicate the fitted curve and the dots 
indicate some values of the normalized pressure amplitude. 
In DSLOW, the dot appears at each third datum. 

D. PROGRAM IMPROVEMENT 

DSLOW was designed to provide three-dimensional graphs. 
For example, the x-axis represents the scaled source 
distance, the y-axis represents the scaled received 
distance, and z-axis represents the normalized pressure 
amplitude. To simplify the presentation, only two- 
dimensional graphs were presented with the x-axis the nor- 
malized pressure amplitude and the y-axis the receiver angle 
6 . All curves are presented with the data fitted with a 
cubic spline. 

The DSLOW program was executed to obtain numerical 
results of the phase angle, the pressure amplitude, and the 
normalized pressure amplitude at each receiver position. 
The first run used double precision for accuracy. 
Difficulties were encountered when the DISSPLA subprogram 
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was attached for making the graphical output, 
precision and DISSPLA were not successful, the 
cision was used, resulting in round-off error. 
3.8 at R 2 = 10.0.) 



When double 
single pre- 
(See Figure 
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(03Q)3TDNV H3M303H 




Figure 3.2 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.3 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.4 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3,5 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.6 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.7 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, R^ varied 
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Figure 3.8 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.9 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.10 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.11 Graphs of receiver angle <$ versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.12 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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Figure 3.13 Graphs of receiver angle 6 versus pressure 
amplitude with R 2 fixed, varied 
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;ed pressure amplitude 



IV. DISCUSSIONS 



A. GRAPHICAL OUTPUT 

The graphs of normalized pressure amplitude as a func- 
tion of receiver angle 6 were investigated for various 
source distances and receiver distances R 2 (Fig. 3.2 to 
3.13) ivhile the other parameters are held contstant. For a 
given wedge angle g and sufficiently small source distance 
(Fig. 3.8 and 3.12), at all receiver distances, the pressure 
increases uniformly towards the bottom. For greater source 
distances, (Fig. 3.9 -3.11) the pressure attains a maximum 
within the wedge for all receiver distance. 

As the receiver distance is increased (Fig. 3.9), a 
pressure minimum develops between the maximum and the 
bottom. An important property of the curves of pressure 
versus receiver angle when there is a maximum and minimum is 
that, at a specific receiver distance, the pressure above 
the minimum can be extrapolated to zero pressure on the 
bottom. (See Fig. 3.10 with R 2 = 32.) This receiver dis- 
tance is called the "transition point." So far, we do not 
know the properties of the transition point. We use the 
transition point for indicating the behavior of the curves 
when the parameter involved is varied. The transition point 
appeared twice in some cases, but in the following 
discussions the first transition point is the only point we 
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will be concerned with. (See Fig. 4.1 for transitions cor- 
respond to R 2 = 4.6 and 6.4.) 

B. GRAPHS CLASSIFICATION 

The development of curves with the source distance (R^) 
and the receiver distance (R 2 ) as variables was observed. 
As R^ or R 2 are varied the curve changes from a linear curve 
to a curve with an observable minimum (Fig. 3.9, R 2 = 5.0) 
and finally to a curve without a minimum (Fig. 3.9, R 2 = 
9.0). Three different types of curves resulted from the 
series of two-dimensional plotting. They are described 
below : 

1 . Type 1 Curves 

Type 1 curves (Fig. 4.2) are those where the sound 
pressure is equal to zero at the surface and maximum at the 
bottom and is almost linearly dependent on depth. These 
curves are most pronounced when the source distance is much 
smaller than the characteristic distance. The closer the 
source is to the characteristic distance, the more nonlinear 
the curves (Figs. 3.3 and 3.4). 

0 

2. Type 2 Curves 

Type 2 curves (Fig. 4.3) are those where the sound 
pressure is zero at the surface, maximum somewhere between 
the surface and the bottom with no minimum. These types of 
curves are generated when the source is placed at a point 
much greater than the characteristic distance. Type 2 
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ANGLE VS. 




Figure 4.1 The plots where there are two transition points, 
R2=4.6 is the first and R2=6.4 is the second 
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curves indicate that the sound energy in the wedge is well 
collimated and that reflection is negligable. 

3 . Type 3 Curves 

Type 3 curves (Fig. 4.4) are those that have a min- 
imum pressure. These curves occur when the source is a dis- 
tance slightly greater than, or less than, the charac- 
teristic distance. Tables 1, 2, and 3 of Appendix D show 
the receiver positions at the first transition points. 
Three different values of B , two different values of 
and two different values of £\/ c 2 were used in making these 
tables. The transition point did not occur when $ = 15°, 
Pl/P 2 = 0.90, c^/c 2 = 1.10. An explanation can be offered 
using the fact that for these particular sound-speed and 
density ratios an angle of intromission exists [Ref. 15]. 
Since the angle of intromission is the grazing angle at 
which the sound energy is completely transmitted into the 
slow bottom, it is plausible that no transition point 
occurs . 

C. TRANSITION POINT 

By varying the wedge angle 8 in small increments A8 = 
0.5° starting with 8 = 5°, and ending at 8 = 7°, it was 

found that transition occurs for source distances within the 
range from 1.0 to 1.5. 

For Ri < 1.0, no transition point was observed; the 
curves are the Type 1. For 1.0 < R^ < 1.5, the evolution of 
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curves as the receiver distance varied can be explained as 
follows: first, the receiver is placed near the source and 
gradually it is shifted further from the source. The mini- 
mum in the pressure decreases reaching the point where the 
curves extrapolate to zero (the first transition point). 
Further detailed observations were made on this particular 
facet by varying the source distance and the receiver 
distances. The results of these observation are tabulated 
arc c raphed in Appendices D and E. When the receiver is 
moved away from the source, the minimum will reach a minimum 
pressure then the pressure increases until it reaches the 
point where the curves again can be extrapolated to zero, 
this is the second transition ooint (See Fig. 4.1). 

For Ri > 1.5, there will be no transition point. The 
curves are the Type 2. 

The transition point as a function of source angle can 
be observed using the tables in Appendix D. In most cases 
the greater the source angle y, the closer the transition 
point is to the apex. Graphs of transition point as a func- 
tion of (Appendix E) indicate that the smaller £5 the more 
regular the curves. This is easy to understand because the 
smaller 0, the more accurate the observation of transition 
point; the greater 0 the less accurate the data. 
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Figure 4.2 Type 2 curves, indicating a pressure 
amplitude nearly linear with depth 
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Figure 4.3 Type 2 curves, indicating a well-collimated 
sound field as the source away from the apex 
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Figure 4.4 Type 3 curves, indicates the presence of 
reflection and refraction near the bottom 
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NORMALIZED PRESSURE AMPLITUDE 



D. PARAMETER VARIATIONS 



Variation of parameters was done by changing one 
parameter of interest while all others were held constant, 
for fixed source and receiver distances, and plotting the 
receiver angle versus normalized pressure amplitude. 

The parameters 3, P 1 /P 2 • an< ^ c l/ c 2' were held constant 
and the pressure amplitude was plotted for various 6, R^ , 
and R 2 • The y's are set at: 3/4, 3/2, and 33/4. Variations 
in the shore distance (Y 0 ) can be made because the program 
is available, but to simplify the investigation, Y 0 was set 
equal to zero for all plots (Fig 4.2 is included as an 
example for Y 0 * 0). 

1. Variations of 3 

Initially, the values of 3 investigated were: 6°, 
10°, and 15°. The najor effect created by altering the 
value of 8 is that, for the sane values of R^ , <?\/ ^ 2 ' an ^ 
c l/ c 2 , the smaller 3, the shorter the transition point (see 
Figs. 3.9 and 3.10, Tables 1, 2, and 3 of Appendix D) . 

2. Variations of y 

The variations of y from y = 3/4 (the source is 
placed near the bottom) to y = 33/4 (the source is placed 
near the surface) are presented in Tables 1, 2, and 3 indi- 
cated that the greater y the shorter the transition point. 
It is not always true, for instance in Table 1 at 3 = 6°, R]_ 
= 1.50, the greater y the longer transition point, for the 
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rises 



( See Appendix D. ) 



3. Variations of c^/c 2 and p^/p 2 

Variations of the acoustical parameters c^/c 2 and 
Pl/P2 were done, but did not give a significant variation of 
the sound pressure profile. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The program DSLOW gives good plots representative of the 
sound energy distribution within the wedge. The sound energy 
can be well collimated by the wedge. This phenomenon is 
strongly affected by the source position. At a source 
position close enough to the apex, sound energy is distri- 
buted linearly with respect to the depth. As the source 
moved away from the apex, the distribution of sound energy 
becomes more complex. Sometimes a minimum is found; this 
minimum may be caused by the presence of sound energy 
reflected by the bottom. 

The source position plays a major role in forming the 
pressure distribution profile. The pressure distribution is 
also very sensitive to the parameter variation at small 
source distances, but it becomes insensitive at large source 
distance. The characteristic distance must have physical 
meanings rather than just an arbitrary number, because when 
the source distance in proximity to the characteristic 
distance, the model is most sensitive. 

The model is restricted when the single precision mode 
generates round-off error and rough curves which do not 
allow for accurate analysis. 
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B. RECOMMENDATIONS 

1. Single precision produces oood results, but failed 
in some cases. Double precision would improve the 
program, but increase the execution tine. This must 
be done by running the program in double precision, 
and accumulating the result in single precision 
before plotting the data by DISSPLA. 

2. It is suggested that the program be run using more 
realistic parameters and observing the effects on 
the characteristic distance and transition point. 

3. Further study validating DSLOW in comparison with 
experimental results is suggested. 
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APPENDIX A 

DSLOW ALGORITHM 



The pressure amplitude calculation 



\| = INT[180,'PJ 


(eqn A.l) 


n 

K ,X = 

1 2tanp tan[arccos(c 2 cj)J 


(eqn A. 2) 


AL = a T ^2 = 0.0001 (constant) 


(eqn A. 3) 


D 2 = Y o 2 + R 1 2 +R 2 2 


(eqn A.4) 


= 2 R j R-> 


(eqn A. 5) 


R 8 = Vd 2 * R 3 c os((N-l)p + y - 5] 


(eqn A. 6 ) 


R 9 =v/D 2 - R 3 cos((N-l)P + y+ 6 ] 


(eqn A. 7) 


S 2 = (.i)INT(Ni' 2 ) 


(eqn A. 8 ) 


W, = 2AL(c 1 ;c 2 ) 2 


(eqn A. 9) 


SI-KRj sin((N-l)p + Y) - 2(INT((N-1):2)P 
+ R 2 sin(2INT((N-l);2)p - 6 ]}|/R g 


(eqn A. 10) 


ci = Vli-si 2 ) 


(eqn A.l 1) 


1=51/0! 


(eqn A.12) 


W 0 = (- c 2 + c i/ c 2 ) 


(eqn A.l 3) 
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(eqn A. 14) 



Y = V\V + W,' 



o 

£ 

II 

N 


(eqn A. 15) 


Y, = V(Y + W 0 ),2 


(eqn A. 16) 


Y 2 = V(Y-W 0 ),2 


(eqn A. 17) 


. _ t,-y 2 

1 (T-Y 2 )- + V 


(eqn A. 18) 


"7 _ Y 1 

2 (T-Y 2 )- + Y,“ 


(eqn A. 19) 


-7 _ (T 2 -y 2 ) 2 - y, 2 

5 (T--Y 2 r T Y/ 


(eqn A. 20) 


2Y 1 T 

6 (T~-Y 2 )“ + y, 2 


(eqn A. 21) 



N 

1 

P, = I (-1) 1NT(N 2) (Z 5 cos(R 8n KjX) + Z 6 sin(R 8n K 1 X}/R 8n (eqn A.22) 
n = 1 



N 

1 

P 2 = I (-1) 1NT(N;2) {-Z 5 sin(R 8n K 1 X) + Z 6 cos(R 8n K,X}/R 8n (eqn A.23) 
n= 1 



N 

1 

P 3 = I (-1) INT(N 2) {Z 5 cos(R 9n K 1 X) + Z 6 sin(R 9n K 1 X}/R 9n (eqn A.24) 

n= 1 
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1 

P 4 = I (*1) INT(NV2) {-Z 5 sin(R 9n K,X) + Z 6 cos(R 9n K,X}/R 9n 
n= 1 


(eqn A. 25) 


p 5 = p, + p 2 


(eqn A. 26) 


P 6 = P 3 + P 4 


(eqn A. 27) 


p tor R iV p 5 J + p 6 2 


(eqn A.28) 
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APPENDIX B 

DSLOW PROGRAM 



************************************************************************ 
** THIS PROGRAM CALLED DSLOW IS CALCULATING THE SOUND PRESSURE ** 
WITHIN THE WEDGE OVERLYING SLOW BOTTOM FLUID AND DOWN-SLOPE ** 



** 



************************** ***Rx^xxxxx* ******* ************ ************ *** 



INTEGER 

REAL’'? 



,S2,N1,J,K,P 



REALM TQQ.T 

pi^=Acosr-i, 

XXK^ 




k k 

kkkkkk 



TggJTQQi.T^, 

7*7<7K7<k'xkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



INPUT .PARAMETERS 

XXXXKXXXXXXXXXXX 



kk 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



B = 

G = 

D = 
Nl = 
R1 = 
“12 = 



RL 
YO = 
D1 = 
CC = 
AL = 
t A = 



WEDGE ANGLE (DEG) 

SOURCE ANGLE TDEG) „ 

RECEIVER ANGLE (DEG) 
b OF IMAGE POINTS 

SOURCE DISTANCE iltj CHARACTERISTIC DISTANCES^ 



S) 



RECEIVER DISTANCE- (IN CHARACTERISTIC DISTANC 
APEX DISTANCE (IN CHARACTERISTIC DISTANCES) 

RHO 1/RHO 2 
C-l/C-2 

* T PH 2 k / Y 0 

-*+*+*-** ■ *t*OF7RECEIVER , POSITION^ 

xXxxxxxxxxxxXxxxxxxxxxxxxxxXxxxxx 

** INITIAL INPUT & RELATIONS 

******X*XXXXX*X*XXX***XXXXX*X*X***************************************** 

c 



kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

kk 



YO = 0.00 
R2 = 9.00 
DO 111 P=l,( 
10.0 
10. Q 
3*B/4 
2* B/4 
l*B/4 
0.0 
0.95 
1.10 
1 .00 
0.75 

i:9J 

1.50 

2.oo 

D.00 

o.oo 

0.30 
1 . 10 
0.0001 



B = 

B = 

G = 

G = 

G = 

D = 

R1 
R1 
R1 
R1 
R2 
R2 
R2 
R2 
R2 
YO 
D1 
CC 
AL_ 

kkkkk ★x kkk i S kkkkkkkkk k k k kk k k k k kkkkkkkkkkkkkk k k k kkk k k k k k k k k k k k k k kkk k kk k k k 

********** 



= 0 

= INT (180 . /B) 

'-M *2 

= TYg*Y^+(Rl*Rl) + (R2*R2) 
=_PI^2 *TAN ( ACOS ( 1 /CC)) *TAN (B ) ) 

u”= t ac6s(i/cc) 



C240 

C 

C 

-241 

-250 

-251 

-270 

-271 

-272 

-800 

-801 



C 

C 



C400 

C 

C 

C 



m miimmihi 

01 = 1/DS0RT( 2 . 0D0 " 



1/DSQ _ 
ORHAT ( ' NOP 

- GE SRR9 
6,240 



* ' I 
WRITE 
FORMAT 
WRITE 
FORMAT 
WRITE 

^FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

* 1 pr 
WRITE 
FORMAT 
WRITE 
DO 10 

D = 



0D00) 



THETA(N) ' ,7X,'IMGE SR R8 ' ,7X, 



B 



»K$>W 



, F5 . 2 , IX , 1 SOURCE ANGLE = 



,F5.2) 

' , F6 . 2 , IX , 



16,241' 

:i.S8i . .... .. 

T SOURCE DISTANCE 5 ' ,F4.2<1X, 'RECEIVER DISTANCE 5 
SHORE DISTANCE 5 ' ,F4.2) 

6,251) R1.R2.Y0 
Y RH01/RH62 5 ' F5.2,5X, 1 C1/C2 5 1 , F5 . 2 , 5X , • ALPHA/K2 5 ' ,F8.4) 

6.270) Dl.Cf 

6.271) T4 



• 2 ) 



6,272 ) 

l£, 'REC^POS^X,.: .RECTANGLE 



iSE ANGLE 

6 , 8 °?y 

6,801) 
M> 0 , A 
M*B/A 



,4#, 'I<I0RM. PRESS' ) 



,2X, 'PRES. AMPLITUDE ' ,3X, 



-110 

-120 



C310 

-20 



DD=B/A 
V =2*B/10 
IF (D .Li .V) 



GOTO 

GOTO 



110 

120 



IF(D.GE.V 
DX 5 D 

D V 5 = DX+(J;1?*B/(10*A) 

I ? ; #°* A > 

d °if?si^gt.o)^ti (N)= (N- l ) *B+G 
IF(Sl.LT.O) T1(N) 5 N*B-G 

51 5 - SI 

R8(N) 5 SORT (D2-R3*COS (T1 (N 
R9 (N) 5 SQRT (D2-R3*COS f T1 (N 
FORMAT (3xll2 ^5X, F5 . 2 . 4X , F6 . 4 ,6 
WRITE (6,310) fl,Ti(N) ,ft8(fl) ,R9(N 
CONTINUE 
PI 5 0.0 
P2 5 0.0 

DO 30 N 5 1,N1 , 

52 5 (-l)**/lNT(N/2)) 

W1 5 2*C2*AL 
II =_INT( (N : l)/2) 




40 i = i.n , , v ^ ^ s 

S(I) 5 AES(R1*SIN(T1(N)-2*I*B) 

rf (^ o ^)+( wlUi )) 
z ) y = z 

Yi '=' qi*$6rt(y+wo) 

Y2 5 -Q1*SQRT(Y-W0) 



Z1 5 T-Y2 
Z2 5 -Yl , 

Z12= CMPLX(Z1R.Z2I) „ 

Z3 5 Zl/Tzl*Zl-l-Z2*Z2) 

Z4 5 -Z2/(Z1*Z1+Z2*Z2) 

Z1 5 t+y2 

Z2 5 Yl 

Z13= CMPLX(Z3R , Z4I ) 

Z5 5 Zl*Z3lZ2*24 
Z6 5 Z1*Z4+Z2*Z3 4 

Z14= CMPLX(Z5R,Z6I) 

E(I) 5 Z5 
F(Il = Z6 

^FORMAT (' W '5X, 1 I ' ,5X, ' ARSIN(S ( I ) *T6 ) 1 ,5X, 'E(I) 1 , 

*WRITE (6 |g3j T ^(l)*E(I) + F(I)*F(I))' ,5X, 'T6*ATAN(F(I)/E(I) ) ' ) 
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C420 

C 

C 

- 40 



FORMAT (3X.I2 < 6X / I2 / 5X.Ft 
^WRITE (6,420) Mil, . AStN I 

CONTINUE 
2 
Z 



X, 12 , 5X. F6 .4 , 10X F6 .4 3X.F 



6.4,3X.F6.4.10X,F7 .4) 
ATAN(F(I)/E(I) ) 



Z S = 



9: 



0 
0 

24 = 0 
Z| = 1 

IF?N,LE:2.0J) 



DO 45 
21 = 
22 = 
Z3 = 

It: 

ze = 

CONTINUE 

Z1 = 
Z2 = 
T = 
Z3 
Z4 



ii) 



E( 

F 

2‘ 

Z6 . 

21*23-22 

zi*24+z; 



Z5 

Z6 

14* 

COS 



GOTO 

II 
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*Z4 

:*Z3 



= zSSla 



N) 



2*24 



C500 


FORMAT ( 


c 


WRITE ( 


C510 


FORMAT ( 


C 


WRITE ( 


C600 


FORMAT ( 


c 




c 


WRITE ( 


- 


DO 60 I 



Z6 = 21*24+22*23 
PI = Pl+S2*Z5/R8 
P2 = P2+S2*26/R8 
' NO, OF I.P ' ,3X 
6,50(5) 

3*. 12 6X,F6 i 4 / 6X £6.4) 



(N) 

!(N) 

f ' O' 
k , IV 



E(REFL)= 1 ,3X, 1 IM(REFL)= ') 



1 1 , SXLliARgll^lf I (f f ) *T6= 1 



, 5X, 'Alto 

' :W‘ 



') 



,3X, 'E(I) 1 , 5X, ’ F ( I ) ' , 5X, 1 EF 1 



I) = SQRT(1-S(I)*S(I)) 



S( 

C('_, , x ,_. 

t = s(i};di vj . , v 

WO = -C2+C(I)*C(I) A 
Y = SgRTUW0*W0) + (Wl*Wl)) 

IF(Y.LE;2)_Y.= 2. 



)RT(Y+WO) 

5Qrt(y-w6 



Q1 K S 
- 01 * 

T-Y2 
-Y~ 

- 22 /( 21 * 21+2 
T+Y2 
Y1 

Z1*Z3-Z2*Z4 
21*24+22*23 



) 



S 0 U ?* 4 * 23 



*1 



2) 



f) ? z 



. 26 
ANl- asin 

)RT 
TAN 



EF = S( 
, AN2= A1 
2X.I2 ‘ 

l§,610J 




I:& 4 4l> F ?(f) 5 ^ F lNt’ 5x ' F6 - 4 - 5x ’ F7 - 4) 
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22 = Z6 . 
T = T4* 
Z3 = COS 
Z4 = -SI 




C700 

C 

C710 

C 

- 30 



:fi° 

-10 



-811 

-31 

C 

-111 






IT (P1*P1+P2*P2 
kN{P2/Plj , x 

' 13 ilX F^l ' 

^ ; i j ^ i i A , r j « ^ , 



p 2 f 6 -“> 



kkkkkk 
kk 
****** 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



FORMAT ( 1 ”LOWE^ PATH HO= '"Sx , 1 RE (REFL)= 1 , 3X , 1 IM(REFL)= ') 

WK 1 1 II l 
FORMAT 
WRITE 
CONTINL. 

K = $+1^ 

DZ(K)=D*T6 

PZ(K)=SQRT(Pl*PltP2*P2)*Rl 
AT(Ki=ATAN(^ / ~" ' 

WRITE (6.8: 

FORMAT (6&, 

CONTINUE 
CONTINUE 
DO 31,L=1,29 
XP=PZ (28 t x , 
pn(l7=pz(l)/xp 
WRITE ( 6, 8ll) L ,D2 
FORMAT (2ft , 13 , 7X , FS 
CONTINUE 
Y0=Y0+0 . 5 
R2=R2*0 . 50 
CONTINUE 
STOP 
END 

rT^:":************************************************************* 

TALL MEDBUF 



PZ (K) , ATAN(P2/P1 } ,PN(K) 

12X, r 7 .4 , IzX, F7 .4 ^ 10X, r7 .4) 



J ! 2 ? 6x* ik ! 6^8 ^h ? .6 !^,F9.6) 



CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 



TEK618 
COMPRS 
NOBRDR 

pageTi 

- AREAZD 
,L HEIGH 




HEIGHT ( . 2 ) 

XNAME( y NORMALIZED PRESSURE AMPLITUDE$ 1 ,29) 
YNAMEf 'RECEIVER ANGLE (DEG) $ ',19 ) 

VTTrvc ( ' ' ' 



YTICKS(5 
XTICKS ( 5 
GRAF (0 . 
GRAF (0 . 
GRAF (0 . 
DOT ; 
GRID(2 
HEADIN, 
MESS AG I 
MESSAGl 
MESSAG 

messag! 

MESSAGl 
MESSAGl 
MESSAGE 
MESSAG 
MESSAG) 
MESSAGl 
MESSAG 
MESSAG) 
MESSAGl 
MESSAGl 
MESSAG 
MESSAGl 
MESSAG) 
MESSAG 
MESSAG 
MESSAG 
MESSAG 
MESSAG 
MESSAG 
MESSAG 
MESSAG 
MESSAG 



6.10.0.50,0. .5.0,15.0) 
10 . ,$ 0 . 0 , 6 . , 6 . 0 , 10 . 0 ) 
3.0 15.0 0. 3.0 15.0) 



i 




.:EC. ANGLE VS .PRESSURES ' ,-100,1.8,1) 
WEDGE Al'IGLEs S' 100, 8.. 7 ) 

15. OS 1 .100, 'ABUf ' , 'ABUf 1 } 

RH01/RA02: $', 100 8.,6.5} x 
~ ~~ •*■100 . 'Abut* , 'Abut’) 

S' .160.8. .6.0? 

106, 'A6 uT ( , 'ABUT') 

ANGlE : S' 100,8. .5, 

,100, ' ABUf ' , 'ABUf r 
100 'ABUT' 'ABUT' 

, 100, 'ABUT' , 'ABUT' 

,100 'ABUT' 'ABUT' 

6lST.: S' 160,8. .5 
,100, 'ABUf' , 'ABUf' ) 

■,100. 'ABUT' . 'ABUT' ) 
iOO , 1 ABUT ' , 'ABUT' ) 

‘IO6, 'ABUT 1 , 'ABUT' 

Sf .INlT: S' .IOO, 8. , 

' ,100, 'ABUT 1 , 'A6 ut‘ 

' ,100 'ABUT' , 'ABUT' 

' ,100, 'ABUT' , 'ABUT' 

'100, 'ABUT' 'ABUT' 

Sf.INtRMT : S', 100, 

',100 . 'abut' , 'Abut 

loo, 'ABUT' , 'ABUT' ) 



7.50 

SOURC 

0.25 

1.10 

4: § s 

§ E Sa D 

1.00 

3.00 

3.00 

REC.D 



>5) 



:> 



L5) 




,4.0) 



' 1UU , nDU 1 , ADU i ) 

C$T . ifilT : S' (IOO. 8 . ,3 
r,100, 'ABUT' , 'A6UT') 



• 5) 
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nnononoonnon 



ST: S' ,100,8. 
MOO, ( ABUf 1 , 



.3.5) 

‘ABUT’ 



CALL MESSAG( 'SHO.D] 

CALL MESSAGf *0.06 
CALL RESET( ^ALL') 

CALL PARA3 ' 

CALL POLY3 
CALL NOCHEK 

CALL CURVE (PN,DZ, 29,0) 

111 CONTINUE V ' ' 

CALL ENDPL(O) 

CALL DONEPL 
STOP 
END 

'k'kit'k'k'k'k'k'k'k'k'k'k'k'kk'k'k-k'k'k+t'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k+i'k'k'k'k'k'k'k-k'kk'k^'k'k'k'k^k.k'k'k'k'k'k'k'k'k-k'k'k 

** For plotting the calculation^ result on TEK618 or Sherpa, erase ** 

** the proper C's in DISSPLA suoproaram above. Put C's in front of ** 

***iii*xS r iIixHI**KlSI*K5x*I^R*i5*5llK*S5S9^5A************************* 
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APPENDIX C 

NUMERICAL RESULTS OF DSLOYV 



WEDGE ANGLE = 10. OC 

SOURCE DISTANCED. 75 
RH01/RH02 3 0.80 
K1X 3 19.44 

REC.POS REC. ANGLE PRES .AMPLITUDE 



SOURCE ANGLE = 2.50 

RECEIVER DISTANCE 3 9.00 
C1/C2 3 1.10 



SHORE DISTANCE 3 0.00 
ALPHA/K2 3 0.0001 



8:?8 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 



10 


0.90 


11 


1.00 


12 


1.10 


13 


1.20 


14 


1.30 


15 


1.40 


16 


1.50 


17 


1.60 


18 


1.70 


19 


1.80 


20 


1.90 


21 


2.00 


22 


3.00 


23 


4.00 


24 


5.00 


25 


6.00 


26 


7.00 


27 


8.00 


23 


9.00 


29 


10.00 



0.007611 
0.007390 
0.007175 
0.006965 
0.006763 
0.006563 
0.006374 
0.006191 
0.006014 
0.005844 
0.00568 
0.00552, 
0.005379 
0.005239 
0.005107 
0.004983 
0.004866 
0.004759 
0.004658 
0.004566 
0.004482 
0.004003 
0.003919 
0.003841 
0.003550 
0.002981 
0.002149 
0.001121 
0.000000 



PHASE ANGLE 

-1.129606 
-1.108071 
-1.085083 
-1.061422 
-1.037255 
-1.011464 
-0.984591 
-0.957120 
-0.928227 
-0.898359 
-0.867409 
-0.835250 
-0.802517 
-0.768184 
-0.732669 
-0.696848 
-0.659586 
-0.622856 
-0.583461 
-0.545011 
-0.505302 
-0.113330 
0.201979 
0.413213 



0.547851 

0.632015 

0.680885 

0.697709 

■0.048596 



NORM. PRESS 



iia 



6.7‘ 

6. Si 

6.403362 
6.215911 
6.03566' 
5 . 85746- 
5.688235 
5.524390 
5.367038 
5.215738 
5.070300 
4.931440 
4.800441 
4.675240 
4.557538 
4.446329 
4.342858 
4.247226 
4.156956 
4.074470 
3.999775 
3.572277 
3.4974" 
3.4277, _ 
3.168519 
2.660608 
1.917873 
1.000000 
0.000029 
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WEDGE ANGLE = 10.00 

SOURCE DISTANCED. 7 5 
RHOl/RHOZ 3 0.80 
K1X = 19.44 



SOURCE ANGLE = 2.50 

RECEIVER DISTANCE 3 4.50 
C1/C2 3 1.10 



REC.POS 



RE C. ANGLE 

Q.OQ 



1.50 

1:^8 

1.80 

1.90 

2.00 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 



SHORE DISTANCE 3 0.00 
ALPHA/K2 3 0.0001 



PRES. AMPLITUDE 

8 : 81 ?lli 

0.036906 
0.036460 
0.036019 
0.035594 
0.035178 
0.034766 
0.034371 
0.03398^ 

0.032867 
0.032511 
0.032163 
0.031823 
0.031492 
0.031166 
0.030846 
0.030532 
0.030226 
0.027368 
0.02461' 
0.0216? 
0.01825 _ 
0.014269 
0.009828 
0.005008 
0.000000 



PHASE ANGLE 

mm 

1.569994 

-1.555159 

-1.538568 

-1.521979 

-1.505386 

-1.488472 

-1.471694 

-1.454889 

-1 .404230 
-1.387417 
-1.370529 
-1.353691 
-1.336905 
-1.320334 
-1.303701 
-1.287225 
-1.270993 
-1.116021 
-0.981865 



-0.872103 
-0.786891 
-0 .72472' 
-0.68275 
-0.660769 
0.273944 






NORM. PRESS 

7.554441 
7 .460239 
7.369551 
7.280349 
7.192361 
7.107576 
7.024406 
6.942234 
' .863304 
.785389 

i iiiiili 

6.562953 
6.491957 
6.422443 
6.3544 ‘ 
6.2833 
6. 2233. 
6.15944 
6 . 09668 
6.035670 
5.464833 
4.917889 
4.325723 
3.644764 
2.84928' 
1.96253. 
1.000000 
0.000001 
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WEDGE ANGLE = 10.00 

SOURCE DISTANCE^. 15 
RH01/RH02 3 0.80 
K1X = 19.44 



SOURCE ANGLE 3 2.50 

RECEIVER DISTANCE 3 2.25 
C1/C2 3 1.10 



SHORE DISTANCE 3 0 . 00 
ALPHA/ K2= 0.0001 



REC.POS REC. ANGLE PRES .AMPLITUDE 



3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 



0.00 
0 . 10 
0.20 
0.30 

8:18 

0.60 

0.70 

0.80 

0.90 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1.90 

2.00 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 



.224816 
. .224336 
0.223838 
0.223305 
0.222733 
0.222113 
0.221472 
0.220805 
0.220096 
0.219349 
0.218573 
0.217751 
0.216893 
0.216000 
0.215070 
0.214095 
0.213090 
0.212046 
0.210963 
0.209839 
0.208671 
0.194767 
0.176765 
0,154712 
0.128893 
0.099744 
0.068024 
0.034481 
0.000000 



PHASE ANGLE 

1.084841 

1.092171 

1.099430 

1.106651 

1.113891 

1.121132 

1.128275 

1.135338 

1.142416 

1.149437 

1.156374 

1.163238 

1.170144 

1.176949 

1.183646 

1.190398 

1.197003 

1.203545 

1.210009 

1.216443 

1.222806 

1.231938 

1.333861 

1.377699 

1.413359 

1.441621 

1.461364 

1.473173 

-1.243458 



NORM. PRESS 

6.519929 

6.505989 

6.491551 

6.476104 

6.459502 

6.441518 

6.422938 

6.403585 

6.383023 

6.361359 

6.338860 

6.315019 

6.290153 

6.264237 

6.237287 

6.209000 

6.179859 

6.149565 

6.118167 

6.085560 

6.051688 

5.648450 

5.126376 

4.486338 

3.738033 

2.892702 

1.972762 

1.000000 

0.000001 
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WEDGE ANGLE = 10.00 

SOURCE DISTANCE=0.75 
RH01/RH02= 0 . T 
KlX = 19. 



SOURCE ANGLE = 2.50 

RECEIVER DISTANCE= 1.13 
Cl/C2= 1.10 



SHORE DISTANCES 0.00 
ALPHA/K2= 0.0001 



REC.POS REC. ANGLE PRES .AMPLITUDE 



10 

11 

12 

13 

if 

16 

17 

18 

19 

20 
21 
22 
23 
2 ' 
2 
2 

2 

28 

29 



0.00 
0.10 
0.20 
0.30 
0.40 
0 . 50 
0 . 60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.30 
1.40 
1.50 
1.60 
1.70 
1.80 
1.90 
2.00 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 



1.900229 

1.923264 

1.944952 

1.965253 

1.984233 

2.001945 

2.018349 

2.033438 

2.047371 

2.060028 

2.071478 

2.081736 

2.090869 

2.098818 

2.105642 

2.111385 

2.115994 

2.119550 

2.122012 

2.123518 

2.123900 

2.074631 

1.939167 

1.730041 

1.459312 

1.138785 

0.780462 

0.396619 

0.000001 



PHASE ANGLE 

-0.125102 

-0.119954 

-0.114742 

-0.109495 

-0.104208 

-o ! sillii 

-0. 082850 
-0.077489 
-0.072146 
-0 . O60/88 
-0.061446 

-8:8l§4ol 

-0.045508 

-0.040234 

-0.034997 

-0.029773 

-0.024576 

-0.019431 

0.029804 

0.074109 

0.112477 

0.144356 

0.169405 

0.187392 

0.198255 

-1.531829 



NORM. PRESS 

4.791066 

4.849142 

4.903825 

4.955010 

5.002865 



, .047521 
5.088881 
5.126925 
5.162055 
5.193968 
5.222836 
5.248699 
5.|71727 

i ! 3089/f 
5.323455 
5.335076 
5.344040 
5.350248 
5.354045 
5.355009 
5.230785 
4.889239 
4.361967 
3.679378 
2.871230 
1.967786 
1.000000 
0.000002 



68 



WEDGE ANGLE = 10.00 

SOURCE DISTANCED.} 5 
RHOI/RHO 27 .O.. 8 O 



K1X = 19.44 



SOURCE ANGLE =2.50 
RECEIVER DISTANCE 3 0.56 
Cl/C2= 1.10 



SHORE DISTANCE= 0 
ALPHA/K2= 0.0001 



REC.POS 

1 



10 

it 

17 

18 

19 

20 
21 



2 

2 . 

26 

27 

28 
29 



REC. ANGLE PRES .AMPLITUDE PHASE ANGLE NORM. PRESS 



0.00 
0 .10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.30 
1.40 
1.50 
1.60 
1.70 
1.80 
1.90 
2.00 
3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 



3.992624 

1:8 if?!? 

S:8!?5§ 7 3 

4.093628 

4.107439 

4.120035 

4.130817 

4.139806 

4.147007 

4.152471 

4-156213 

4 . 158244 

4.158599 

4-157314 

4.154387 

4.149846 

4.143713 

4.136028 

4.126792 

3.953794 

3.647243 

3.22468 

2.70326, 

2.100711 

1.435872 

0.728641 

0.000000 



0 . 3064^ 



. .31466. 
0.318822 
0.323001 
0 .327198 
0.331284 
0.335505 
0.339736 
0.343969 
0.348202 

8 : 11111 ! 

8:38%! 

0.369249 

0.373416 

0.377565 

0-381693 

0.335796 

0.389876 

0.428945 

0.464082 

0.494513 

0.519774 

0.539597 

0.553331 

0.562396 

-0.847966 



5.479548 
5.512648 
5.543023 
5.570723 
5.595763 
5.618168 
5.637122 
5.654409 
5.669208 
5.681543 
5.691426 
5.698924 
5.704061 
5.706848 
5.707335 
5.705572 
5.701555 
5.695323 
5.636906 
5.676359 
5.663683 
5 . 426256 
5.005542 
4.425613 
3.710009 
2.883053 
1.970616 
1.000000 
0.000000 



.00 
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WEDGE ANGLE 3 10.00 

SOURCE DISTANCED. 75 
RH01/RH02 3 0.80 
K1X 3 19.44 



SOURCE ANGLE = 2.50 

RECEIVER DISTANCE 3 0.28 
C1/C2 3 1.10 



SHORE DISTANCE 3 0.00 
ALPHA/K2 3 0.0001 



REC.POS REC. ANGLE PRES .AMPLITUDE 



8 

10 

li 

13 

11 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 



0.0 

0 . 1 . 

0.20 

0.30 

0.40 

0.50 

8:58 
0.80 
0.90 
1.00 
1 . 10 
1.20 
1.30 
1.40 
1.50 
1.60 
1.70 
1.80 
1.90 
2.00 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 



0.642356 

0.636151 

0.629941 

0.623727 

0.617507 

0.611286 

0.605059 

0.598325 

0.592588 

0.586344 

0.580099 

0.573846 

0.567591 

0.561331 

0.555062 

0.548791 

0.542518 

0.536235 

0.529949 

0.523657 

0.517360 

0.454118 

0.390374 

0.326155 

0.261508 

0.196416 

8 * JiieSI 

0.000001 



PHASE ANGLE 

-0.270501 
-0.270035 
-0.269567 
-0.269111 
-0.268658 
-0.268203 
-0.267760 
-0.267313 
-0.266874 
-0.266441 
-0.266014 
-0.265585 
-0.265165 
-0.264745 
' 264339 
. 263923 
-0.263521 
-0.263120 
-0.262721 
-0.259016 
-0.255795 
-0.253081 
-0.250935 
-0.247807 
-0.246445 
-0.245626 
-0.026186 



NORM. PRESS 






9.60207. 
9.507351 
9.412549 
9.31771' 

1 ! 12^778 
9.032702 
8.937539 
8.842341 
8.747032 
8.651687 
8.556263 
8.460711 
8.365122 
8.269497 
8.173734 
8.077916 
7.982002 
7.886027 
6.922033 
5.950404 
4.971522 
3.986116 
2.99393' 
1.99845. 
1.000000 
0.000009 
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APPENDIX D 
TABLES 

TABLE 1 

RECEIVER DISTANCE AT THE FIRST TRANSITION POINT, 
FOR CONSTANT p, p 2 = 0.S0, Cj c 2 = 1.10 



P = 6°, K 1 X= 32.61 





R, = 1.10 


R, = 1.20 


o 

m 

II 

tC 


R, = 1.40 


R,- 1.50 


II 

"CD 


5.2 


5.0 


5.9 


7.3 


9.8 


Y = P 2 


5.2 


4.7 


5.7 


6.8 


9.9 


Y=3P4 


4.0 


4.5 


5.4 


7.0 


10.5 



P= 10°’ K ,X = 19.44 





R, =0.80 


R, = 0.90 


R, = 1.00 


R, = 1.10 


R, = 1.20 


Y = P 4 


17.5 


24.0 


33.0 


52.0 


72.0 


Y = P 2 


12.5 


17.0 


24.0 


42.0 


60.0 


Y= 3p'4 


10.6 


17.0 


22.0 


40.0 


58.0 



P= 15°’ K ,X= 12.79 





R, = 1.30 


'w' 

II 


R, = 1.50 


R, = 1.60 


R, = 1.70 


II 

•CD 


8.1 


8.2 


8.7 


9.8 


11.4 


Y = P 2 


6.9 


7.0 


7.8 


8.9 


10.9 


II 

Ci 

"CD 


9.8 


12.0 


15.0 


19.0 


29.0 



71 



TABLE 2 



RECEIVER DISTANCE AT THE FIRST TRANSITION POINT, 
FOR CONSTANT p t p : = 0.80, Cj.c,- 1.20 

P = 6°, K,X= 22.53 





R, = 1.30 


R, = 1.40 


R, = 1.50 


R, = 1.60 


R, = 1.70 


ll 

-C3 

-U 


3.37 


3.6 


3.92 


4.4 


4.96 


y = P 2 


3.155 


3.42 


3.78 


4.36 


4.96 


y=3P 4 


3.07 


3.34 


3.74 


4.2 


5.1 



P = 10°’ K,X= 13.43 





R, = 0.80 


R, =0.90 


R, = 1.00 


R, = 1.10 


R, = 1.20 


y = p 4 


6.1 


11.3 


18.0 


40.0 


60.0 


II 

"CD 


6.5 


8.2 


14.5 


25.0 


40.0 


y=3P4 


5.05 


5.45 


12.0 


18.0 


30.0 



P= 15°’ K |X = 8.84 





R, =0.80 


R| = 0.90 


R, = 1.00 


R, — 1.10 


o 

II 

& 


II 

"CD 


36.0 


40.0 


52.0 


60.0 


64.0 


II 


24.0 


26.0 


32.0 


46.0 


54.0 


Y=3p 4 


19.50 


23.0 


28.0 


38.0 


58.0 
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TABLE 3 



RECEIVER DISTANCE AT THE FIRST TRANSITION POINT, 
FOR CONSTANT p t p : = 0.90. Cj c 2 = 1.10 



p= 6°’ K,X = 32.61 





R, = 1.10 


R, = 1.20 


R, = 1.30 


R, = 1.40 


R, = 1.50 


y = p4 


3.66 


4.1 


4.9 


6.2 


11.0 


II 

*C3 


3.35 


3.9 


4.7 


7.0 


no 


y=3p'4 


3.24 


3.76 


4.75 


7.3 


9.S 



P = 10°, K,X= 19.44 





R, = 1.30 


R, = 1.40 


R, = 1.50 


R, = 1.60 


R| = 1.70 


Y = P 4 


9.8 


11.0 


14.0 


40.0 


60.0 


y=P 2 


7.6 


8.6 


11.0 


19.0 


40.0 


Y = 3P4 


6.7 


7.6 


10.0 


21.0 


45.0 
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R1 VERSUS R2 AT TRANSITION POINT 



APPENDIX E 



GRAPHS OF Rj VERSUS R 2 AT THE FIRST TRANSITION POINT 




Figure E.l Rj vs R 2 at the first trans. points, for P = 6°, Pj/P 2 = 0.80, C|/C 2 — l.IO. 
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B=6(DEG),RH0 l/RH02=.8,CJ/C2= 1 . 1 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.2 Rj vs R 2 at the first trans. points, for P= 10°, P 1 /P 2 = 0.80, Cj/C 2 = 1.10. 
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B=10(DEG),RH01/RH02=.8,Cl/C2=l.l 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.3 Rj vs R 2 at the first trans. points, for P= 15°, Pi/p 2 -0.80, Cj/c 2 — 1-10. 
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B=15(DEG),RH01/RH02=.8,C1/C2=1.1 




R1 VERSUS R2 AT TRANSITION POINT 




Figure E.4 Rj vs R 2 at the First trans. points, for P= 6°, pj/p 2 = 0.80, Cj/cj- 1-20. 
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B=6(DEG).RH01/RH02=.9.Cl/C2=l.l 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.5 Rj vs R-> at the first trans. points, for P= 10°, P j/P2 — 0*80, c j,^- l 



78 



B=10(DEG),RH01/RH02=.9,C1/C2=1.1 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.6 Rj vs R 2 at the first trans. points, for P= 15°, P 1 /P 2 = 0.80, c j / c 2 = 1-20. 
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B=15(DEG),RHO1/RHO2=.0,C1/C2=1.2 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.7 Rj vs Ro at the first trans. points, for P= 6°, P 1^2“ 0.90, C|/c 2 - 1.10. 
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B=6(DEG),RH0 1/RH02=.8,C l/C2=l .2 



R1 VERSUS R2 AT TRANSITION POINT 




Figure E.8 Rj vs R 2 at the first trans. points, for p= 10°, P j /P 2 ~ 0.90, Cj/C 2 = 1.10. 
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B=1O(DEG),RHO1/RHO2=.0,C1/C2=1.2 
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